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Abstract 

In  this  paper,  we  propose  a  model  for  dynamic  magne¬ 
tostrictive  hysteresis  in  a  thin  rod  actuator.  We  derive 
two  equations  that  represent  magnetic  and  mechanical 
dynamic  equilibrium.  Our  model  results  from  an  appli¬ 
cation  of  the  energy  balance  principle.  It  is  a  dynamic 
model  as  it  accounts  for  inertial  effects  and  mechani¬ 
cal  dissipation  as  the  actuator  deforms,  and  also  eddy- 
current  losses  in  the  ferromagnetic  material. 

We  also  show  rigorously  that  the  model  admits  a  pe¬ 
riodic  solution  that  is  asymptotically  stable  when  a  pe¬ 
riodic  forcing  function  is  applied. 

1  Introduction 

There  is  growing  interest  in  the  design  and  control  of 
smart  structures  -  systems  with  embedded  sensors  and 
actuators  that  provide  enhanced  ability  to  program  a 
desired  response  from  a  system.  Applications  of  interest 
include:  (a)  smart  helicopter  rotors  with  actuated  flaps 
that  alter  the  aerodynamic  and  vibrational  properties  of 
the  rotor  in  conjunction  with  evolving  flight  conditions 
and  aerodynamic  loads;  (b)  smart  fixed  wings  with  ac¬ 
tuators  that  alter  airfoil  shape  to  accomodate  changing 
drag/lift  conditions;  (c)  smart  machine  tools  with  actua¬ 
tors  to  compensate  for  structural  vibrations  under  vary¬ 
ing  loads.  In  these  and  other  examples,  key  technolo¬ 
gies  include  actuators  based  on  materials  that  respond 
to  changing  electric,  magnetic,  and  thermal  fields  via 
piezoelectric,  magnetostrictive  and  thermo-elasto-plastic 
interactions. 

Typically  such  materials  exhibit  complex  nonlinear 
and  hysteretic  responses  (see  Figure  1  for  an  example  of 
a  magnetostrictive  material  Terfenol-D  used  in  a  com¬ 
mercial  actuator).  Controlling  such  materials  is  thus  a 
challenge.  The  present  work  is  concerned  with  the  de¬ 
velopment  of  a  physics-based  model  for  magnetostrictive 
material  that  captures  hysteretic  phenomena  and  can  be 
subject  to  rigorous  mathematical  analysis  towards  con¬ 
trol  design. 

In  Section  2  we  propose  a  model  for  magnetostric¬ 
tion  that  describes  the  dynamic  behaviour  of  a  thin  rod 
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tures  URI  Contract  No.  DAAL03-92-G0121  and  under  the  MURI97 
Program  Grant  No.  DAAG55-97-1-0114  to  the  Center  for  Dynamics 
and  Control  of  Smart  Structures  (through  Harvard  University). 


actuator.  It  is  a  low  (6)  dimensional  model  with  10  pa¬ 
rameters  and  hence  is  suitable  for  real-time  control.  The 
model  incorporates  features  observed  in  a  commercial  ac¬ 
tuator  [1],  like  the  hysteretic  behaviour  of  magnetostric¬ 
tion  as  a  function  of  the  external  field;  dependence  on 
the  rate  of  the  input,  eddy  current  losses,  inertial  effects 
and  mechanical  damping  effects.  In  Section  3  we  ana¬ 
lyze  this  model  for  periodic  forcing  functions.  Using  the 
Schauder  Fixed  Point  Theorem,  we  prove  that  the  solu¬ 
tion  is  an  asymptotically  stable  periodic  orbit,  when  the 
parameters  are  subject  to  certain  constraints. 

2  Thin  magnetostrictive  actuator  model 

We  are  interested  in  developing  a  low  dimensional  model 
for  a  magnetostrictive  actuator.  The  main  motivation  is 
to  use  it  for  control  purposes.  Therefore  the  starting 
point  of  our  work  is  Jiles  and  Atherton’s  macroscopic 
model  for  hysteresis  in  a  ferromagnetic  rod  [2].  In  our 
model,  we  treat  the  actuator  itself  along  with  the  asso¬ 
ciated  prestress,  magnetic  path,  to  be  a  mass-spring  sys¬ 
tem  with  magneto-elastic  coupling.  As  we  show  later, 
our  model  is  only  technically  valid  when  the  input  signal 
is  periodic.  However,  this  is  the  case  in  many  applica¬ 
tions  where  one  obtains  rectified  linear  or  rotary  motion 
by  applying  a  periodic  input  at  a  high  frequency  to  these 
actuators.  For  instance  in  our  hybrid  motor  [3,  4],  we 
produced  a  rotary  motion  using  both  piezoelectric  and 
magnetostrictive  motors  in  a  mechanical  clamp  and  push 
arrangement. 

In  an  earlier  work,  we  explored  the  connections  be¬ 
tween  a  bulk  ferromagnetic  hysteresis  model  and  energy 
balance  principles  [5].  We  present  in  this  paper,  an  ex¬ 
tension  of  this  theory  to  include  magnetostriction  and 
eddy  current  losses.  This  is  done  by  equating  the  work 
done  by  external  sources  (both  magnetic  and  mechani¬ 
cal)  ,  with  the  change  in  the  internal  energy  of  the  mate¬ 
rial,  change  in  kinetic  energy,  and  losses  in  the  magneti¬ 
zation  process  and  the  mechanical  deformation. 

SWbat  +  6Wmech  =  6Wmag  +  SWmagd  +  SWd 

Change  in  internal  energy 

+  SLmag  +  SLel+  SK  (1) 

V - V - ' 

losses  Change  in  kinetic  energy 

In  Equation  1,  SK  is  the  work  done  in  changing  the 


kinetic  energy  of  the  system  consisting  of  the  magnetoe- 
lastic  rod,  SWmag  is  the  change  in  the  magnetic  potential 
energy,  SWmagei  is  the  change  in  the  magnetoelastic  en¬ 
ergy,  SWei  is  the  change  in  the  elastic  energy,  SLmag  are 
the  losses  due  to  the  change  in  the  magnetization,  and 
5Le[  are  the  losses  due  to  the  elastic  deformation  of  the 
rod. 

The  expression  for  SWmag  will  be  given  shortly.  The 
elastic  energy  is  given  by  We|  =  ^Ix1 ,  where  x  is 
the  total  strain  multiplied  the  length  of  the  actuator. 
Chikazumi  [6]  derives  an  expression  for  the  magnetoe¬ 
lastic  energy  density  of  a  three  dimensional  crystal.  It 
is  of  the  form  where  the  strain  components  multiply  the 
square  of  the  magnetization  components.  In  our  one 
dimensional  case,  we  can  similarly  write  down  the  fol¬ 
lowing  expression  for  the  magnetoelastic  energy  Wmagei  ■ 

Wmagel=bM'2xV 

where  b  is  the  magneto  elastic  coupling  constant  and  V 
is  the  volume  of  the  magnetostrictive  rod.  M  is  the  av¬ 
erage  magnetic  moment  of  the  rod.  The  expression  for 
the  magnetic  hysteresis  losses  SLmag  is  due  to  Jiles  and 
Atherton.  The  motivation  for  this  term  is  the  observa¬ 
tion  that  the  hysteresis  losses  are  due  to  irreversible  do¬ 
main  wall  motions  in  a  ferromagnetic  solid.  They  arise 
from  various  defects  in  the  solids  and  are  discussed  in 
detail  by  Jiles  and  Atherton  [2] .  The  mathematical  con¬ 
sequences  of  this  hypothesis  is  discussed  in  detail  in  our 
earlier  work  [5]. 


SL 


mag 


V  k  sign(H)  (1  —  c)  dMirr 


The  line  integral  implies  that  the  integration  is  carried 
out  over  one  full  cycle  of  the  input  voltage/current  which 
is  assumed  to  be  periodic.  The  reason  for  this  will  be 
discussed  shortly.  The  losses  due  to  mechanical  damping 
are  assumed  to  be  SLei  =  <f  c\  x  dx.  The  change  in  the 
kinetic  energy  SK  =  j>  meff  x  dx.  Therefore, 


6Wbat 


/j.qH  d.M 


Let  an  external  force  F  in  the  x  direction  produce  a 
uniform  compressive  stress  in  the  x  direction  a  within 
the  actuator.  The  total  displacement  of  the  edge  of  the 
actuator  rod  be  x.  Thus  the  mechanical  work  done  by 
the  external  force  in  a  cvcle  of  magnetization  is  given  by 

[7], 

=  F  dx 

The  total  work  done  by  the  battery  and  the  external 
force  is 


6Wbat  +  SWmech  =  V 


[it)  II  d.M  T 


F  dx 


We  see  that  adding  the  integral  of  any  perfect  differential 
over  a  cycle  does  not  change  the  value  on  the  left  hand 
side.  Therefore, 


8Wbat+8Wmech  =  V  (J>  p.0H  d.M  + 
Equations  2  and  3  give, 


aM  d\l  + 


F  dx 


(3) 


V  no  §  {H  +  aM  dM  ••• 

+  f  (F  —  lx  —  Clx  —  meffX  —  VbM2)  dx  =  SWmag 

+V  j>  ksign(H)  (1  —  c)  <IM-, (4) 


Define  the  effective  field  to  be, 


He 


=  H  +  aM  — 


2  b  Mx 
Ik) 


As  the  integration  is  over  one  cycle  of  magnetization,  we 
have 


He  dM  =  - 


M  d.He 


“I”  fyiech  —  “I”  y  TFleff  %  dx  .  .  . 


6K 


+  V  (p  b M2  dx  +  V  (p  2b M x  dM  +  (b  lx  dx 


SWmagel  6Wet 

+  V  <b  k  sign(H)  (1  —  c)  d.Mjrr  +  j)  a  x  dx  (2) 


SL„ 


SLe 


It  was  observed  in  [5] ,  that  if  M  is  a  function  of  He  then 
there  are  no  losses  in  one  cycle.  This  is  the  situation  for 
a  paramagnetic  material  where  M  =  Man  is  given  by 
Langevin’s  expression  as  a  function  of  He.  Hence  for  the 
lossless  case,  the  magnetic  potential  energy  is  given  by, 

SWmag  =  ~V  j  Man  d.He 
Thus  Equation  4  can  be  rewritten  as 


Now  we  obtain  expressions  for  the  left  hand  side  of  the 
above  equation.  For  a  thin  cylindrical  magnetostrictive 
actuator,  with  an  average  magnetic  moment  M,  and  an 
uniform  magnetic  field  in  the  x  direction  H ,  the  work 
done  by  the  battery  in  changing  the  magnetization  per 
unit  volume,  in  one  cycle,  is  given  by 


V/Jo  $  ( Man  -  M  -  k^n(H)(i-0  d^)  dHe 

+  j>  (F  —  d.  x  —  c i  x  —  mef  f  x  —  b  M 2  V)  dx  =  0 

Note  that  the  above  equation  is  valid  only  if  //.  M.  x,  x 
are  periodic  functions  of  time.  In  other  words,  the  tra¬ 
jectory  is  a  periodic  orbit.  We  now  make  the  hypothesis 


that  the  following  equation  is  valid  when  we  go  from  one 
point  to  another  point  on  this  periodic  orbit. 


V/jo  f  (Man  -M-  dMj^)  dHe 

+  J  (F  —  dx  —  ci  x  —  meff  x  —  b  M2  V)  dx  =  0 

The  above  equation  is  assumed  to  hold  only  on  the  peri¬ 
odic  orbit.  Since  dx  and  dHe  are  independent  variations 
arising  from  independent  control  of  the  external  prestress 
and  applied  magnetic  field  respectively,  the  integrands 
must  be  equal  to  zero. 


Fed 


N 2  87 TP  B'2  -  A2 
l  m  B2  +  A2 


The  actual  work  done  by  the  battery  in  changing  the 
magnetization  and  to  replenish  the  losses  due  to  the  eddy 
currents  in  one  cycle  is  now  given  by 


“1“  Peddy  dt 

=  -V  I IqM  dHe  +  Peddy  dt 


Ma 


k  sign(H)  (1  —  c)  dMirr 
fo  dHe 

mef  f  x  +  ci  x  +  dx  +  b  M2  V 


0  (5) 

F  (6) 


Figure  2  shows  a  schematic  of  the  full  model.  The 
hysteretic  inductor  stands  for  the  magnetostrictive  actu¬ 
ator  model. 


Jiles  and  Atherton  relate  the  irreversible  and  the  re¬ 
versible  magnetizations  as  follows  [2] , 


Qualitative  analysis  of  the  magne 
tostrictive  actuator  model 


A  f  —  l\  fre  v  Mjrr. 

Ad  rev  —  C  ( Afan  Adjrr ). 


rev 

dM  dMirr  dMan 

——  =  6m  (I  -  C)  -rrr-  +  c- 


dH 

where  8m  is  defined  by. 


dH 


dH 


(7) 


8m 


0  :  H  <  0  and  Man(He)  -  M(H)  >  0 
0  :  H  >  0  and  Man(He )  -  M(H)  <  0 
1  :  otherwise. 


(8) 

Finally  after  some  algebraic  manipulations,  the  equa¬ 
tions  for  the  magnetostriction  model  are  given  by 


dM 

dt 


^  -AM  dt  ^ 

mef  f  x  +  ci  x  +  dx  +  b  M2  V  =  F  (10) 


It  is  very  important  to  note  that  the  model  equations  (9 
-  10)  are  only  valid  when  all  the  M,  H  are  periodic  in 
time.  What  we  mean,  is  that  solution  of  the  equations 
represent  the  physics  of  the  system  under  these  condi¬ 
tions.  But,  usually  in  practice  we  do  not  know  apriori 
what  state  the  system  is  in.  Then  can  we  use  the  above 
model?  The  answer  in  the  affirmative  is  provided  in  this 
section.  We  show  analytically  that  even  if  we  start  at 
the  origin  in  the  M  —  H  plane  (which  is  usually  not  on 
the  hysteresis  loop),  and  apply  a  periodic  input  H ,  we 
tend  asymptotically  towards  this  periodic  solution. 

It  was  shown  in  an  earlier  work  [5],  that  Equation  (9) 
has  an  orbitally  asymptotically  stable  limit  cycle  when 
6  =  0  (no  coupling)  and  the  input  is  co-sinusoidal.  The 
situation  is  much  more  complicated  when  the  coupling  is 
non-zero.  It  remains  to  be  shown  that  there  exists  an  or¬ 
bitally  asymptotically  stable  limit  cycle  for  co-sinusoidal 
inputs  H,  with  6  yf  0. 


The  inputs  to  the  above  set  of  equations  are  an(j 
F,  while  x  is  the  mechanical  displacement. 

A  magnetostrictive  material  has  finite  resistivity,  and 
therefore  there  are  eddy  currents  circulating  within  the 
rod.  Using  Maxwell’s  equations,  we  can  derive  the  fol¬ 
lowing  simple  expression  for  the  power  losses  due  to  eddy 
currents  [3]. 

V2  lm  b2+a 2 
eddv  N2  87 tp  B 2  -  A2 

where  A ,  B  are  the  inner  and  the  outer  radii  of  the  rod, 
lm  is  its  length,  A7  is  the  number  of  turns  of  coil  on  the 
rod,  and  V  is  the  voltage  across  the  coil  of  the  induc¬ 
tor.  Hence  the  eddy  current  losses  can  be  represented 
equivalently  as  a  resistor  in  parallel  with  the  hysteretic 
inductor.  This  idea  is  quite  well  known  and  a  discussion 
can  be  found  in  [6]  or  [3] .  From  the  above  expression  for 
the  power  lost,  the  value  of  the  resistor  is, 


3.1  The  uncoupled  model  with  periodic  pertur¬ 
bation 

Before  studying  the  full  coupled  system,  we  consider  the 
effect  of  periodic  perturbations  on  the  uncoupled  models. 
Define  state  variables, 

x  i  =  H 

X2  =  M 
Vi  =  x 
y-2  =  x 

Let, 

•Cl  +  - x'2 

a 

Then  the  state  equations  are: 


Theorem  2  Consider  the  system  given  by  Equations 
x1  =  u.  (11)  (11-14),  with  input  given  by  Equation  (16)  and  6^0. 

If  (aji,a;2)(0)  =  (0,0),  then  the  0,-limit  set  of  the  system 
is  an  asymptotically  orbitally  stable  periodic  orbit. 

x-2  =  f2{xi,X2,x3,x4,g{t))u  (12) 


where  the  function  /2(-)  is  obtained  by  substituting  the 
state  variables  in  Equation  (9)  and  g(-)  for  x(-). 


X3  =  sign(u). 


(13) 


X4  = 


X3  <  0  and  coth(^)  —  \  —  ff-  >  0 


X3  >  0  and,  cotli(  %)  —  A  —  jf-  <  0 
otherwise. 

(14) 


y  =  Ay  +  h2(t) 


0 


1 


(15) 


1 

where  y  =  [  yi  y2  ]  ;  A  = 

g{ •),  h(-)  are  Ajl  periodic  functions.  The  input  is  given 
by, 


u(t)  =  U  cos  {cot).  (16) 


3.1.1  Analysis  of  the  uncoupled  magnetic  sys¬ 
tem 

The  proof  of  existence  and  uniqueness  of  trajectories  for 
the  system  (11  -  14)  is  exactly  as  in  the  earlier  paper  [5], 
with  the  some  modifications. 

Theorem  1  Consider  the  system  of  equations  (11  -  If) 
with  b  0.  Let  the  input  be  given  by  Equation  (16). 
Suppose 


\g\  <  G  (17) 

and  a  =  a  —  ALU  safisfles 

MO  J 


Mo 


(1 


c  a  Ms 
3  a 

a  c  Ms  \ 
3  a  ) 


<  1 

-2  d  Mg  >  0. 


(18) 

(19) 


Also 


0  <  c  <  1.  (20) 

Then  the  exists  a  solution  to  the  system  with  initial 
condition  x(0)  =  (0,0).  Moreover  this  solution  is  unique 
for  all  time  t  >  0  and  lies  in  the  compact  set  [— w,  w]  x 

Proof  The  proof  is  very  similar  to  that  of  the  system 
with  6  =  0  [5]. 

□  . 

From  now  on  until  the  end  of  this  section,  it  is  always 
assumed  that  the  parameters  satisfy  conditions  (17-20). 


Proof 

The  proof  is  identical  to  the  one  with  6  =  0  [5] . 

□ 


Denote  the  periodic  solution  of  the  perturbed  mag¬ 
netic  system  (11  -  14)  with  perturbation  g(-)  and  in¬ 
put  u(-),  as  x(-).  It  is  a  two  dimensional  vector  and 
a  T  =  —  periodic  function.  Define  the  sets  B  = 
{</>  G  C([0,T\,n)  :  H  <  ft;  \4>(t)  -  < 

M1\t-t\  V  t,t  e  [0, T]},  D  =  {ip  e  C{[0,T],1Z)  : 
H\  <  Ih;  -  if{t)\  <  M-2  \t  -  f\  V  t,f  G  [0, T]}, 
where  $% ,  fd-2 ,  Mi,  M2  are  positive  constants.  Let  Vi ,  V2  ■ 
C([0,T],7^2)  — y  C([0,T],7^)  denote  the  projection  oper¬ 
ators  defined  by  Vi  (f,g)  =  f  and  V2  (f,g)  =  g. 

Consider  the  mappings  Q  :  B  — ^  C([0, T], g(- )  h-> 
x(-)  and  hi:  D  ^  £([0,T],7?.2);  h(-)  ^  y(-).  We  first 
show  Q  to  be  continuous. 

Theorem  3  Q  is  a  continuous  map. 

Proof  Let  the  system  (11  -  14)  be  represented  by 

x  =  f(t,x,a)  ;  (t,  x)  £  D  C  TZ 3 

where  a  =  a  —  AklUl  anc[  jj  js  an  open  set.  The  state 
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x  is  2-dimensional  because  the  discrete  states  X3  and  X4 
are  functions  of  x\ .  x-2  and  u.  Let  the  initial  condition 
be  (xi,x2)(0)  =  (0,0). 

If  ff n  —>■  g  in  the  uniform  norm  over  [0,T]  where  T  is 
the  period  of  /,  then  an  — >  a.  Consider  the  sequence  of 
systems  x  =  fn{t,x)  =  f{t,x,an).  As  /  is  continuous  in 
a,  fn  — >  /  in  the  uniform  norm  if  dn  — >  a  (Theorem 
8).  The  solutions  of  each  of  the  systems  {/„}  and  / 
exist  and  is  unique  for  t  £  [0,T].  Then  by  Theorem  9, 
the  solutions  (pn(t)  of  x  =  fn{t,x)  converge  uniformly  to 
(j>{t)  the  solution  of  x  =  f(t,,x,a)  for  t,  £  [0,T], 

Consider  the  time  interval  [T,  2T].  We  have  shown  that 
<j>n(T)  — >  <p{T).  Then  again  by  Theorem  9,  <pn(t)  — »• 
<j>{t)  for  t  £  [T,2T],  Thus  we  can  keep  extending  the 
solutions  <f>„  (t)  and  <j>(t)  and  obtain  uniform  convergence 
over  any  interval  [ mT ,  (to+1)T]  where  m  >  0.  Therefore, 
for  each  m  and  e  >  0,  there  exists  N(m)  >  0  such  that 
|^n  -  <p\  <  |  V  n  >  N(m). 

By  Theorem  2  there  exist  asymptotically  orbitally  stable 
periodic  orbits  xn  of  the  systems  x  =  fn(t,x)  and  x  of 
the  system  x  =  f(t,x,ct).  Hence  for  each  e  >  0,  there 
exists  M  >  0  such  that  \xn  —  <pn\  <  |  and  \x  —  <p\  < 
|  V  m  >  M  and  t,  £  [mT,(m+l)T], 

Hence  for  all  n  >  N(M )  and  t  £  [mT,  (m  +  1)T]  where 
m  >  M,  we  have  \xn  —  x\  <  \xn  —  <pn\  +  \<p„  —  <p |  + 
\x  —  4>\  <  e.  Hence  Q  is  a  continuous  map. 

□ 


3.1.2  Analysis  of  the  uncoupled  mechanical  sys¬ 
tem 

In  this  subsection,  we  consider  the  mechanical  system 
with  periodic  perturbation  given  by  Equation  (15).  We 
assume  the  homogenous  system  (that  is,  (15)  with  h(t)  = 
0)  to  be  asymptotically  stable.  The  relevant  results  are 
collected  in  the  appendix. 

Theorem  4  Consider  the  system  (15).  If  the  eigenval¬ 
ues  of  A  have  negative  real  parts  and  h(-)  is  an  peri¬ 
odic  function,  then  (15)  has  an  Jj-  periodic  solution  that 
is  asymptotically  07'bitally  stable. 

Proof  This  follows  from  Lemma  10  and  Theorem  11  in 
the  appendix. 

□ 

Theorem  5  If  the  eigenvalues  of  A  have  negative  real 
parts,  then  B  is  a  continuous  map. 

Proof  This  again  follows  from  Lemma  10  and  Theorem 

11. 

□ 

3.1.3  Analysis  of  the  coupled  magnetostriction 
model 

In  this  section,  we  prove  the  existence  of  an  orbitally 
asymptotically  stable  periodic  orbit  for  the  magne¬ 
tostriction  model. 

Let  D\  denote  the  range  of  P2  °  Q  and  B i  denote  the 
range  of  V\  o  B.  Thus  V-2  °  Q  :  B  h->-  Di  and  V\  o  B  : 
D  n-  Bi. 

Theorem  6  There  exists  ab  >  0  such  that  if  |6|  <  b 
then  P2  o  Q  '■  B\  D\  and  V\  o  B  :  D\  n-  B\ . 

Proof  First  we  show  that  the  sets  B\  and  D\  have  the 
same  structure  as  that  of  B  and  D  respectively.  Then 
we  choose  b  so  that  the  domains  and  ranges  of  Q  and  B 
are  suitably  adjusted.  Choose  /3i  =  Ms  and  Al\  =  U 
in  the  definition  of  the  set  D. 

By  Theorem  1,  the  elements  of  Di  are  uniformly 
bounded  by  Ms.  Let  x  =  Q  g.  Therefore  V2  °  Q  g  =  X2- 
Now  ,r2(t2)  - S2(ti)  =  J0  X2{ti  +s  (to  — ti)) (t2 -t\)  ds 
by  the  Mean  Value  Theorem.  As  the  parameters  of  the 
system  (11  -  14)  satisfy  the  conditions  (17  -  20),  the 
vector  field  f(t,x)  u(t )  is  uniformly  bounded.  Therefore 
\x-i (ti)  —  rc2(t2) |  <  Mi  | to  —  t\\.  Thus  D\  has  the  same 
structure  of  D. 

Let  y  =  Bh.  Therefore  yi  =  V\  o  B  h.  The  elements 
of  B 1  are  uniformly  bounded  because  B  is  linear  in  /r 
and  the  functions  h  G  D  are  uniformly  bounded.  \y\ \  < 
\y\  <  (Pi  o  B\  M|  =  02-  We  need  to  choose  b  so  that 
a  =  a  —  AMi  defined  in  Theorem  1  satisfies  Conditions 
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(18)  and  (19).  Such  a  non-zero  b  obviously  exists.  Now 
yi(h)~yi(h)  =  Jo  yi(ti+s(t2-ti))(t2-ti)  ds  by  the 
Mean  Value  Theorem.  \y\  <  \A\fii  +  bV 0l  =  Mi- 
Therefore  |yi(t2)  —  yi{ti)\  <  Mi  |i2  —  ti|.  Thus  f?i  has 
the  same  structure  of  B. 


Our  choice  of  b  >  0  ensures  that  if  |6|  <  b  then 
V2  o  Q  :  Bi  Di  and  P\  o  B\  D\  i->  Bi . 

□ 

We  now  return  to  the  dynamic  model  of  magne¬ 
tostriction  (9,10)  and  prove  the  main  theorem  of  this 
paper. 

Theorem  7  Consider  the  dynamic  model  for  magne¬ 
tostriction  given  by  Equations  (11  -  16).  Suppose  the 
matrix  A  has  eigenvalues  with  negative  real  parts  arid 
the  parameters  satisfy  conditions  (18-20)  with  the  mag¬ 
netostriction  constant  b  <  b  defined  in  the  statement  of 
Theorem  6.  Then  there  exists  an  orbitally  asymptotically 
stable  periodic  orbit  of  the  system. 

Proof  The  sets  B 1  and  Di  are  compact  and  convex  by 
Theorem  12.  Then  B 1  x  Di  is  compact  in  the  uniform 
product  norm  by  Theorem  13.  Obviously  it  is  also  con¬ 
vex. 

Let  be  defined  as,  $  :  Bj  x  Di  — >  Bi  x 
Di\^  (x2,yi)  =  {Vi  o  B(x2),P2  °Q(yi))-  Then  ^  is 
continuous  because  P2  o  Q  and  Vi  o  B  are  continuous  by 
Theorems  3  and  5,  and  the  continuity  of  the  projection 
operator. 

Then  by  the  Schauder  Fixed  Point  Theorem  (Theo¬ 
rem  14),  there  exists  a  limit  point  of  the  mapping  in 
the  set  Bi  x  Di.  This  gives  us  the  periodicity  of  the  two 
state  variables  ;c2  and  y  1 .  In  general,  the  fixed  point  may 
not  be  unique,  but  when  the  initial  state  is  the  origin,  the 
O  limit  set  is  unique  by  the  uniqueness  of  solutions.  Now, 
(2/1, 2/2)  =  Qx 2  and  by  Theorem  4,  (j/i,j/2)  is  an  asymp¬ 
totically  stable  periodic  orbit.  Also  (.ti,:e2)  =  Byi  and 
by  Theorem  2,  (xi,x2)  is  an  asymptotically  stable  peri¬ 
odic  orbit.  The  other  state  variables  (^3,^4)  are  periodic 
because  they  are  determined  by  Xi,  xo  and  u. 

□ 

A  Mathematical  Preliminaries 

Theorem  8  If  X  and  Y  are  normed  linear  spaces  and 
f  is  a  mapping  from  X  to  Y  ,  then  f  is  continuous  at  x 
if  and  only  if  for  each  sequence  {xn}  in  X  converging  to 
x  we  have  {f(xn)}  converging  to  f(x)  in  Y . 

Theorem  9  Suppose  {/„},  n  =  1,  2,  •  •  • ,  is  a  sequence 
of  uniformly  bomided  functions  defined  and  satisfying  the 
Caratheodory  conditions  on  ail  open  set  D  in  1Zn+1  with 
linjn-Hx.  fn  =  fo  uniformly  on  compact  subsets  of  D. 
Suppose  ( tn,xn )  is  a  sequence  of  points  in  D  converging 
to  ( to,xo )  in  D  as  n  — >  00  and  let  <pn(t),  n  =  1,  2,  •  •  •, 
be  a  solution  of  the  equation  x  =  fn(t. ,  x)  passing  through 
the  point  (tn,  x„).  If  is  defined  on  [ a,b ]  and 

is  unique,  then  there  is  an  integer  no  such  that  each 
<j>n(t ),  n  >0,  can  be  defined  on  [a,b]  and  converges  uni¬ 
formly  to  4>o(t)  uniformly  on  [a,b]. 

Consider  the  homogenous  linear  periodic  system 
x  =  A(t)x 


(21) 


and  the  non-homogenous  system 

x  =  A(t)x  +  f(t)  (22) 

where  A(t  +  T)  =  A(t),  T  >  0  and  A(t)  is  a  continuous 
n  x  n  real  or  complex  matrix  function  of  t. 

Definition  1  If  A(t)  is  annxn  continuous  matrix  func¬ 
tion  on  (  —  00,00)  and  V  is  a  given  class  of  functions 
which  contains  the  zero  function,  the  homogenous  sys¬ 
tem  x  =  A(t)x  is  said  to  be  noncritical  with  respect  to  V 
if  the  only  solution  of  Equation  (21)  which  belongs  to  V 
is  the  solution  x  =  0.  Otherwise,  system  (21)  is  said  to 
be  critical  with  respect  to  V. 

The  set  'Pi  denoting  the  set  of  T-periodic  continuous 
functions  is  a  Banach  space  with  the  sup-norm.  That  is, 
I/I  =  sup-00<t<00  \f(t)\;  f  G  VT-  Let  B  denote  the  set 
of  continuous  bounded  functions  from  1Z  to  lZn. 

Theorem  10  [8]  (a)  System  (21)  with  Aft )  G  Vt  is 
noncritical  with  respect  to  B  if  and  only  if  the  character¬ 
istic  exponents  of  (21)  have  nonzero  real  parts. 

(b)  System  (21)  with  A  G  Vt  is  noncritical  with 
respect  to  Vt  if  and  only  if  I  —  X (T)  is  nonsingular, 
when  X(t ),  .Y(0)  =  I,  is  a  fundamental  matrix  solution 

of  (21). 

Theorem  11  Suppose  A  is  in  Vt-  Then  the  nonho- 
mogenous  equation  (22)  has  a  solution  ICf  in  Vt,  if  and 
only  if  system  (21)  is  noncritical  with  respect  to  Vt-  Fur¬ 
thermore,  if  system  (21)  is  noncritical  with  respect  to 
Vt,  then  ICf  is  the  only  solution  of  (22)  in  Vt  and  is 
linear  and  continuous  in  f. 

Theorem  12  Suppose  D  is  a  compact  subset  of  lZm; 
M,  3  are  positive  constants  and  A  is  the  subset  of 
C(D,lZn)  such  that  <f>  G  A  implies  \<j>\  <  (3 ;  | <p(t)  — 

<)>(t)\  <  M  1 1  —  t\  for  t,t  G  D.  Then  the  set  A  is  convex 
aiid  compact. 

Theorem  13  [9]  Let  A,  B  be  compact  subsets  of  X .  •  || ) 
aiid  ( Y ,  |  •  1 2 )  respectively.  Then  A  x  B  is  compact  (under 
either  of  the  standard  metrics  || (x,y),  {x,y)\\  =  \x—x\\  + 
\y-y\-2,  \\(x,y),(x,y)\\-2  =max.(\x  -  x\i,\y -  y\i)).) 

Theorem  14  (Schauder)  [8]  If  A  is  a  convex,  compact 
subset  of  a  Banach  space  X  and  f  :  A  —t  A  is  continu¬ 
ous,  then  f  has  a  fixed  point  in  A. 
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Figure  1:  ETREMA  MP  50/6  Actuator  displacement  (Microns)  vs  current  (Amps)  characteristic  at  different  driving 
frequencies  . 


Figure  2:  A  thin  magnetostrictive  actuator  in  a  resistive 
circuit. 


